Shear- Viscosity to Entropy Density Ratio of a Relativistic Hadron Gas 
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Ultrarelativistic heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) are thought to 
have produced a state of matter called the Quark-Gluon-Plasma, characterized by a very small shear 
viscosity to entropy density ratio ry/s, near the lower bound predicted for that quantity by Anti- 
deSitter space/Conformal Field Theory (AdS/CFT) methods. As the produced matter expands and 
cools, it evolves through a phase described by a hadron gas with rapidly increasing r)/s. We calculate 
ri/s as a function of temperature in this phase and find that its value poses a challenge for viscous 
relativistic hydrodynamics, which requires small values of rj/s throughout the entire evolution of 
the reaction in order to successfully describe the collective flow observables at RHIC. We show that 
the inclusion of non-unit fugacities will reduce rj/s in the hadronic phase, yet not sufficiently to be 
compatible with viscous hydrodynamics. We therefore conclude that the origin of the low viscosity 
matter at RHIC must be in the partonic phase of the reaction. 



Ultrarelativistic heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) are thought to have produced 
a Quark Gluon Plasma (QGP) with the characteristics of 
a near ideal fluid[i|, S, ll 3- 

One of the most important 
current challenges in QGP research is to quantify the 
transport coefficients of this novel state of matter. Re- 
cently, attention in the field has been primarily focused 
on the shear viscosity to entropy density ratio ry/s. Cer- 
tain supersymmetric gauge theories with gravity duals Q 
suggest a lower bound of 77,„i,i = s/Att for this quantity, 
often referred to as the KSS bound [6|. Relativistic vis- 
cous hydrodynamic calculations require very low values 
of 77/s in order to reproduce the RHIC eUiptic flow (^2) 
data 0, B 0|- However, current calculations assume a 
fixed value of rj/s throughout the entire evolution of the 
system and neglect its temperature dependence. The ex- 
act value of r]/s in these calculations is only known within 
a factor of w 3, due to systematic uncertainties related 
to the choice of equation of state and initial conditions 
used [1, [1^ . A viscous hydrodynamical analysis [q] finds 
that i]/s should lie within a range of 0.08-0.24 depending 
upon the choice of initial conditions and equation of state. 
This finding is supported by calculations of r]/s for pure 
gluonic QCD, which yield values close to the KSS bound, 
and indirect estimates of rj/s from calculations of the dif- 
fusion of heavy quarks, elliptic fiow measurements, and 
transverse momentum correlations that arrive at roughly 
comparable values 11, 13, 13, 14. 15|. 



Note that the shear viscosity of matter in a relativistic 
heavy ion collision is a time- dependent quantity. While 
the partonic phase of such a collision is expected to have 
a very low value of 77/5, after hadronization occurs rj/s 
is expected to rapidly increase. In order to quantita- 
tively constrain the viscosity of the deconfined phase of 
a relativistic heavy ion collision, a separate calculation 
of the hadronic viscosity is necessary. Several investiga- 
tions suggest rj/s should reach a minimum in the vicinity 
of a phase transition or crossover [3,[l3|- In particular, 
[13] argue that ij/s should decrease as a function of ris- 



ing temperature in the hadronic phase, and then increase 
in the deconfined phase. They suggest that rj/s should 
reach a minimum near the deconfinement transition, but 
also remark that the perturbative methods used in their 
arguments are not applicable near Tc- It has also re- 
cently been argued that the existence of Hagedorn states 
[ist will strongly decrease hadronic 77/ s, especially close 
to Tc- It should be noted though that all calculations 
performed thus far have assumed kinetic and chemical 
equilibrium, equivalent to unity light quark fugacities. 
Such an assumption may be reasonable for the formation 
of a QGP in a relativistic heavy ion collision at RHIC; yet 
the different timescales of chemical and kinetic freeze-out 
in the hadronic phase of the reaction imply an acquiring 
of non-unit particle species dependent fugacities as the 
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20|,l2l|,l22] 



system evolves in the hadronic phase 

Several analytic calculations of 77 and rj/s for sim- 

t hadronic systems have previously been performed 
H m. Hi, 113, Hi]. These analytic calculations 
solved the linearized Boltzmann equation, in which the 
cross sections in the collision integral were treated using 
different techniques, such as chiral perturbation theory, 
effective NN theory, and phenomenological amplitudes. 
However, even the most sophisticated analytic calcula- 
tions include one or two hadronic species at most, and 
such a binary mixture clearly is a very crude approxima- 
tion of the hadronic matter present at RHIC. Sophisti- 
cated Monte Carlo microscopic transport models, which 
include the full range of the hadronic spectrum (including 
resonances) [29], provide a far more realistic description 
of the hadronic matter created in relativistic heavy ion 
collisions. Several studies within those models have fo- 
cused on equilibration and thermodynamic properties of 
infinite hadronic matter [SOJ, and have extracted trans- 
port coefficients of hadronic gases [3l|, H^, albeit none 
have performed a systematic study of rj/s. 

In this Letter, we use a microscopic transport model 
known as the Ultrarelativistic Quantum Molecular Dy- 
namics (UrQMD) model, described in [iol . Hi], to sim- 
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FIG. 1: The viscous correlator for T=67.9 MeV. The relax- 
ation time associated with this plot is 136 fm/c. 



ulate infinite equilibriated hadronic matter. We confine 
the particles comprising the system to a box with peri- 
odic boundary conditions in coordinate space [30| . and 
the cohisions force the system into equilibrium. We ver- 
ify that the system has achieved chemical equilibrium by 
checking whether the particle multiplicities in our sys- 
tem saturate as a function of time, and comparing such 
yields to an independent statistical model of a hadron res- 
onance gas (SHARE) ^] . We verify kinetic equilibrium 
by checking the momentum distributions of the hadrons 
in our system for isotropy, and fitting particle spectra to 
Boltzmann distributions. 

In order to extract the shear viscosity of our system, 
we employ the Kubo formalism. In addition, we compute 
the entropy of the system and evaluate ry/s as a function 
of temperature and baryo-chemical potential. The Kubo 
formalism relates linear transport coefficients to near- 
equilibrium correlations of dissipative fluxes and treats 
dissipative fluxes as perturbations to local thermal cqui- 



librium '3a 
viscosity is 

V = 



The Green-Kubo formula for shear 



£r / dt{7T-y{6,0)7r-y{r,t)),^^n, (1) 



where T is the temperature of the system, t is the post- 
equilibriation time (the above formula defines i = as 
the time the system equilibriates) , and tt^^ is the shear 
component of the energy momentum tensor tt'"'. The 
expression for the energy momentum tensor tt^" is 



.M- - / d'p^fix,p), 



(2) 



where f{x,p) is the phase space density of the particles 
in the system. Our system does not assume any inter- 
particle potential, and UrQMD treats the hadrons as 
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FIG. 2: In the left panel, pressure versus energy density. The 
speed of sound extracted from the slope yields = 0.18. In 
the right panel, entropy density versus temperature. 



point particles uniformly distributed in coordinate space, 
which implies 



TT " = — 



N 

E 



(3) 



where V is the volume of the system. In the Kubo for- 
mula, the averaging symbol denotes an averaging over the 
ensemble of events generated in our simulation. A repre- 
sentative sample of the correlations of the shear compo- 
nent of the energy momentum tensor is given in Fig. 1. 
The correlation functions are empirically found to decay 
exponentially in time, hence we assume an exponential 
ansatz to integrate the correlation function over time. 

In order to compute the entropy of the system, we use 
the Gibbs formula s = (^ ^+p-^bpb ^ _ r^]^^ energy density, 

pressure, baryo-chemical potential, and number densities 
of the relevant chemical species are extracted once the 
system has equilibriated. The pressure is computed via 



P ^ 



3V ^ 



bf (j) 
P°(j) 



(4) 



and the chemical potentials are extracted by calculat- 
ing ratios of the particle yields. Extracting the entropy 
accurately from a microscopic transport model requires 
some thought; the medium cannot be simply treated as 
an ideal gas of massless particles, and one should note 
that the specific entropy contributions from particles of 
different masses are not the same. We verify our cal- 
culation for s via a scaling relation between the entropy 
density and the temperature of an equilibriated system 

at fixed volume: s T''= , with Cg the speed of sound. 
The left frame of Fig. 2 shows pressure as a function 
of energy density. The speed of sound can be extracted 
from = (^), evaluated along an adiabat; the slope 
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FIG. 3: The viscosity to entropy ratio in chemical and kinetic 
equilibrium. Also shown is the result of rj/s for chiral pions, 
and a pQGP of 3 massless quarks. 



from the left frame in Fig. 2 yields — 0.18. This value 
has been fitted to the extracted entropy densities as a 
function of temperature and hence satisfies the entropy 
scaling relation. In addition to calculating the entropy 
via the Gibbs formula, we have calculated the entropy by 
determining the individual particle species multiplicities 
in the system and summing over them, weighted with 
their specific entropies (such specific entropies are ex- 
tracted via the aforementioned statistical model). Since 
both methods for the entropy extraction agree to within 
5-10%, and the entropy scaling relation is verified, we 
conclude that our entropy computation is accurate and 
represents a hadronic medium including multiple differ- 
ent particle species. 

Using the aforementioned techniques, we present 77/ s 
as a function of temperature in full equilibrium for zero 
baryo-chemical potential in Fig. 3. Also illustrated in 
that figure is the calculation of rj/s for chiral pions 25 1 
and 3 flavor perturbative QCD |38[. Our results are in 
qualitative agreement with the result for chiral pions - 
they suggest that the minimum value of rj/s in a, hadron 
gas should occur near Tc [17\ . The minimum value found 
for rj/s for the equilibrium and zero baryo-chemical po- 
tential case is ~ 0.9, significantly higher than the KSS 
bound of rj/s « 0.08. If the minimum value of ry/s for 
hadronic matter in the range of hadronic freezeout in- 
deed occurs at that value, our results of rj/s would pose 
a serious problem for the application of viscous hydro- 
dynamics to the hadronic collision evolution at RHIC, 
since a value of rj/s of at most 0.24 is needed to repro- 
duce RHIC elliptic flow data [8] . 

However, calculating rj/s in full kinetic and chemical 
equilibrium, as has so far been common practice, may be 
unsatisfactory: while a statistical model analysis of parti- 
cle yields and ratios at RHIC indicate a chemical freeze- 
out temperature in the vicinity of Tchem ~ 160 MeV 
(391 . [4Q] . hydrodynamic calculations indicate that a ki- 



netic freeze-out temperature of Tkin ~ 130 MeV is re- 
quired to describe the momentum distributions of final 
state hadrons [13, 21 , [H, |4lj . The separation of time- 
scales (and hence temperatures) for chemical and kinetic 
freeze-out imply that the hadronic phase of a relativis- 
tic heavy ion reaction actually acquires increasingly non- 
unit fugacities as the system cools after chemical freeze- 
out and evolves out of eauilibriuni [2 0|.| 2 ll|. One should 
note that microscopic calculations [l9|,l22| of the hadronic 
evolution do not require explicit introduction of such fu- 
gacities, since the respective phenomena are the result 
of dynamically changing collision rates for inelastic vs. 
elastic hadronic rescattering processes. In the context 
of our hadron gas calculations, a non-unit fugacity can 
be induced by initializing the system with a surplus of 
pions (kaons) relative to the \Tr,K = exp(/j,^^i<-/r) — 1 
value for the corresponding temperature and performing 
the viscosity measurement before the system relaxes into 
chemical equilibrium. A similar effect can be obtained 
by inducing a finite net-baryon density, which results in 
a finite baryon-chemical potential. 

Figure 4 displays our calculation for rj/s as a function 
of temperature at non-unit fugacities or at finite baryo- 
chemical potential. The largest effect is seen for finite 
baryochemical potentials. The reduction in r]/s can be 
understood classically if one associates a non-unit fugac- 
ity with an increasing particle density of that correspond- 
ing species in the system. Since ?7 ~ § (with p the mean 
momentum of the particle), increasing the particle densi- 
ties of a given species in the system will lead to a reduced 
mean free path, which in turn reduces the viscosity (in 
the case of adding baryons, the average cross section will 
rise, with similar effect). Similarly, increasing the mul- 
tiplicities of different species in the system will enhance 
the entropy density, hence rj/s will decrease. Taking non- 
unit fugacities into account, our values for rj/s in the 
range of hadronic chemical freezeout can be as low as 
Tj/ s ~ 0.4 — 0.5. However, this is still significantly above 
the suggested values rj/s ~ 0.08 — 0.24 from viscous hy- 
drodynamics calculations, and we should keep in mind 
that the hot hadronic matter created at RHIC, though 
not at unit pion or kaon fugacities, still is most likely at 
As « 1. This observation enables us to constrain the 
origin of the low viscosity phase in a relativistic heavy 
ion reaction. While a perturbative QGP is expected to 
have a large value oi ij/s, it should strongly decrease as 
a function of decreasing temperature in the strongly cou- 
pled non-perturbative region, and its minimum should 
be reached at or near Tc on the deconfined side of Tc- 
•q/s may then exhibit a sharp rise or even a disconti- 
nuity as it crosses Tc from T ~ T+. Discontinuities in 
rj/s at Tc have been computed in weakly-coupled scalar 
field theories, and arguments have been given to suggest 
the discontinuity in 77/ s at a phase transition could be a 
universal feature of a much larger class of systems [l6j . 

In summary, we have calculated the viscosity over en- 
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FIG. 4: The shear viscosity and viscosity to entropy ratio as 
a function of temperature for different values of ^ and Xt^^k- 
Note the strong decrease in - for large values of the fugacities. 



tropy density ratio rj/s of a hadron gas as a function 
of temperature, baryo-chemical potential and fugacities. 
We have demonstrated that the inclusion of non-unit par- 
ticle fugacities, which are bound to arise due to the sepa- 
ration of chemical and kinetic freeze-out during the heavy 
ion collision evolution, will reduce the value of rj/s, but 
not to the value necessary to ensure the successful ap- 
plication of viscous hydrodynamics to the full collision 
evolution at RHIC. Our calculation of 77/s in a hadron 
gas from a microscopic transport model therefore con- 
strains the origin of the low viscosity matter produced in 
a relativistic heavy ion collision, which must occur in the 
deconfined phase, possibly in the range 1 < T/Tc < 2. 
At the formation of the hadronic phase, which is thought 
to occur in chemical equilibrium, rj/s will experience a 
sharp increase. However, subsequently its value may de- 
crease again due to the system evolving out of chemical 
equilibrium. Near kinetic freeze-out tj/s will rise with de- 
creasing temperature. This lends credence to the notion 
that the dynamics of the evolution of a collision at RHIC 
is dominated by the deconfined phase exhibiting very low 
values of rj/s. 
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